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^ Abstract. The dust formed in extended circumstellar envelopes of long-period vari- 

O ables and Miras has a strong influence on the envelope dynamics. A radiatively driven 

instability caused by the formation of dust leads to the development of an autonomous 
dynamics characterised by a set of distinct frequencies. We study the interplay between 
the envelope's internal dynamics and an external excitation by a pulsating star. 

4^ 1. Introduction 

Q Long-period variables (LPVs) and Miras are radially pulsating, highly evolved stars 

on the asymptotic giant branch (AGB). Their cool, extended atmospheres are excellent 
^ sites for the formation of complex molecules and dust particles. The interplay between 

■ dust formation and stellar radiation results in circumstellar envelopes (CSEs) generat- 

ing slow mass loss, which finally enriches the interstellar medium with processed mate- 
^ rial. The observed light curves exhibit irregularities over intervals of several pulsation 

QQ cycles, which are at least partially caused by the phenomenon of dust formation. 

To investigate the complex dynamical behaviour of carbon-rich envelopes in detail, 
^ CSEs can be considered as nonlinear multi-oscillatory systems, whose eigenfrequencies 

and normal modes are controlled by the intrinsic timescales of various coupled physical 
and chemical processes. This is done by applying established methods of non-linear 
dynamics such as Fourier Analysis and the study of stroboscopic and Poincare maps 
upon the results of self-consistent model calculations ( [Fleischer et al. |1992| ). 
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^ 2. Eigendynamics 

^ The envelopes around highly luminous AGB stars, e.g. = 3 x 10"^ L©, - \ M©, 

= 2450 K, and CjO - 1.25 (henceforth model A) develop a self-maintaining os- 
cillation caused by dust formation (exterior /c-mechanism) even without the additional 



input of mechanical momentum by an underlying stellar pulsation (cf. [Fleischer et al. 
1995[|Hofner et al. |1995| ). In their power spectr um (Fig.[T} left) one can clearly see the 
eigenfrequency ^ 2450 d"^ and its overtones ( Dreyer eFal. [2009 )). 



In the case of lower stellar luminosities, e.g. = 9 x 10^ L©, - \ M©, 
= 2500 K, and CjO - 1.75 (henceforth model B), the system generates no self- 
induced oscillation. An additional input of energy and momentum is needed to deter- 
mine the CSE eigendynamics. We excite such an envelope by white noise. Thereby 



we adopt a normal noise distribution with zero mean and intensity a (see Dreyer et al. 
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2010[ ). The resulting power spectrum (Fig.[TJ right) exhibits different eigenmodes. The 
position of the minima (marked in grey) coincides with the dust nucleation zone located 
at r ^ 2.25 R^. This confirms that the dust is responsible for creating the characteristic 
frequencies of the envelope. 





Figure 1 . Power spectra of the radial outflow velocity u for a high luminous LPV 
purely caused by the exterior A:-mechanism (model A, left) and caused by a stochastic 
excitation for a standard luminous LPV (model B, right). The position of the dust 
formation zone in the envelope is highlighted. 



3. Interaction with external excitation 

Once the eigenfrequencies have been found, the question arises as to how these modes 
interact with an external excitation caused by the stellar pulsation. For this purpose, 
we studied the response of our reference models for a series of harmonic excitations 
with different periods P and strength Aw. Figure [2] shows the most dominant envelope 
frequencies over the excitation period at constant strength. For model A (left) we can 
identify three different dynamic domains, namely eigenmode -dominated, irregular and 
pulsation-dominated domain. The same study for model B (right) reveals a similar 
behaviour. In contrast to model A we have not found an eigenmode dominated domain. 

Figure [3] shows the maps of our reference models for a typical oscillation period of 
LPVs. In the maps it can be seen that for a constant phase angle (black) the system tends 
to stay inside a finite number of stripes. Similar to the orbital resonance in celestial 
mechanics, we can find integers / and j which fulfil the resonance condition iP = jPcsE- 
This suggest that the systems returns to the same physical state after / excitation periods 
P or j resonance periods Pcse- For model A, we find that / = 27 and j = 4. For model 
B the CSE response is double periodic, which means that the envelope period equals 
two times the excitation period ( [Dreyer et al. |2009l|2010[ ). 



4. Synthetic Hght curves 

Figure |4] presents synthetic light curves for model A obtained for various excitations at 
several wavelengths. The light curves in the first and the last panels show monoperi- 
odic behaviour but with differing periods. For the unexcited CDS (left) the dominant 
timescale is the dust formation scale P^ = 2450 d, whereas in the case of the highly- 
excited model (right), this timescale is determined by the stellar pulsation P = 1500 d. 
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Figure 2. Most dominant frequencies of a CSE /cse normalised to its eigenfre- 
quency (triangles) excited with various pulsation periods P near the dust nucle- 
ation zone for model A (left) and model B (right). The set of CSE-eigenmodes is 
depicted by circles, the eigenperiod is labelled with P^. The excitation frequencies 
(solid line) and harmonics (dash-dotted lines) are also shown. 




Figure 3. Maps of the (u,p,T) phase space in the dust nucleation zone for 
model A (left) and model B (right). The stroboscopic maps were obtained 
by sampling [u(n0.02P),p(n0.02P),T(n0.02P)] (grey) and the Poincare maps by 
[u(nP),p(nP), T(nP)] (black) for 1 < ^ < LW/^J,^ e N. 



In the curves in the second panels the eigenperiod of the CSE P^ = 2450 d and the exci- 
tation period P = 388 d interfere with each other. Nevertheless, the dominant timescale 
is the slightly shifted dust formation timescale Pcse ~ 2700 d. The light curves in the 
third panels show double-periodic behaviour. The eigenperiod is detuned to the first 
harmonic of the excitation period Pcds = 2P = 3000 d. 



5. Special application: IRC +10216 



Winters et al. | ( |1996| ) have found that the fundamental parameters: = 24 x 10^ L©, 
= 2500 K, C/0 = 1.2 as well as P = 650 d, and Au = 6kms-^ well 

So, we use this set of parameter to model 
To determine the envelope 



M, = 1 Mo, n 

reproduce observations of IRC +10216. 
and investigate the dynamical behaviour of the envelope, 
eigendynamics the investigation follows the approach used for CSEs around standard 
and low luminous AGB stars. Since the central star pulsates with a well defined pe- 



riod of P = 650 d (e.g. |Menten et al. ||2006| ), we excite the model exactly with this 
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Figure 4. Synthetic light curves of model A for various excitations at different 
wavelength. The magnitudes are related to the mean magnitude. 



frequency. As can be seen in the power spectrum (Fig. |5| left) the CSE is domi- 
nated by the pulsation frequency and its first harmonic close to the star. With onset 
of dust formation at r 2.5 these frequencies disappear and the envelope frequency 
/cDS = 0.26 X 10"^"^ ^ (3900 d)"^ controlled by the timescale of the dust forma- 
tion process gains strength. The corresponding Poincare map (Fig. |5} right) shows four 
clearly distinct stripes (black). A resonance such as in the case of the previously studied 
models A and B is not found. 




Figure 5 . Power spectrum of velocity u (left) and maps of the (u, p, r)-phase space 
(right) of the periodically excited IRC -hi 02 16 model. The excitation frequency is 
highlighted in the spectrum. 



Figure|6]compares specific observations and synthetic light curves in different NIR 
bands. The stellar pulsation period P = 650 d can be seen in all bands. An additional 
period of six times the stellar pulsation period is also visible (dashed-dotted line). This 
is the dust-determined eigenperiod of the CSE Pcse ~ 3900 d. Due to the short obser- 
vation period the data cannot reproduce the long envelope period. Nevertheless on can 
see a continuous upward trend in the observation data (J, H, K-band). This provides a 
hint of the underlying dust-determined period. 
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Figure 6. Comparison between calculated (line) & observed (points) light curves 
of IRC +10216 at different NIR- wavelengths. The magnitudes are related to the 
mean magnitude. The CSE-eigenperiod Pqse is marked with a dashed-dotted line. 



6. Summary 

We investigated the dynamics of carbon-rich envelopes around AGB stars. We showed 
that the presence of dust has an important influence on the dynamical behaviour of 
the entire system. Firstly, it amplifies the momentum coupling between matter and 
radiation. Secondly, it enforces an additional dynamical behaviour (exterior k- mech- 
anism) to the dust envelope. The consequence of the coupling is the development of 
an autonomous envelope dynamics which is in fact excited by the stellar pulsation but 
caused by the timescales of physical processes in the envelope. This dynamics can be 
characterised by distinct frequencies which do not need to be identical with the stel- 
lar pulsation frequency. Depending on the excitation the CSEs show distinct dynamics 
regimes. We have found a resonance condition which relates the stellar pulsation pe- 
riod to the envelope modes. Besides the pulsation period, light curves exhibit timescales 
which are determined by the dust. 
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